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Abstract 

Many theories of the early universe predict the existence of a multiverse where bubbles continu¬ 
ously nucleate giving rise to observers in their interior. In this paper, we point out that topological 
defects of several dimensionalities will also be produced in de Sitter like regions of the multiverse. 
In particular, defects could be spontaneously nucleated in our parent vacuum. We study the evo¬ 
lution of these defects as they collide with and propagate inside of our bubble. We estimate the 
present distribution of defects in the observable part of the universe. The expected number of such 
nearby defects turns out to be quite small, even for the highest nucleation rate. We also study 
collisions of strings and domain walls with our bubble in our past light cone. We obtain simulated 
full-sky maps of the loci of such collisions, and find their angular size distribution. Similarly to 
what happens in the case of bubble collisions, the prospect of detecting any collisions of our bubble 
with ambient defects is greatly enhanced in the case where the cosmological constant of our parent 
vacuum is much higher than the vacuum energy density during inflation in our bubble. 
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I. INTRODUCTION 


A wide variety of particle physics models predict the existence of a multiverse, where 
bubbles of different vacua nucleate and expand in an eternally inflating background. Occa¬ 
sionally bubbles collide, and observers inside the bubbles can search for the imprints left by 
these collisions on the cosmic background radiation (CMB). Detecting such imprints would 
be tremendously important, as they would provide a direct evidence for the existence of 
a multiverse. The dynamics of these collisions have been studied analytically as well as 
numerically over the years M and the analysis of their observational signatures has now 
become an active area of research; see, e.g., mn] and references therein. 

In the present paper, we shall discuss another potentially observable manifestation of a 
multiverse. It has been shown in na that topological defects, like spherical domain walls and 
circular loops of cosmic string can be spontaneously produced in a de Sitter-like universe. 
The initial radii of walls and strings are close to the de Sitter horizon. Once they are 
formed, they are stretched by the inflationary expansion, resulting in a scale-invariant size 
distribution of walls and/or strings. If indeed we live in an expanding bubble, our bubble 
will collide with topological defects nucleating in the inflating parent vacuum. Any such 
collisions that occurred within our past light cone can give rise to potentially observable 
effects. 

The physical properties of the defects depend on the embedding vacuum and may change 
signihcantly as the defects cross from the parent vacuum to the bubble interior. In fact, 
matching defect solutions may not even exist inside the bubble. Collisions of defects with 
our bubble would then result in the disintegration of the defects. The opposite limit is 
when the helds making up the defects do not interact with the held of the bubble, except 
gravitationally. Then the defects crossing from parent vacuum into the bubble will only 
respond to the change of spacetime geometry. Walls and strings encountered by our bubble 
will then remain, at least partially, in the bubble interior and may be observable directly. 
These defects, originating in the false vacuum, can in principle be much heavier than strings 
and walls that can be formed in phase transitions inside our bubble. 

Our goal in this paper is to study the dynamics of the intruding defects and to estimate the 
number of potentially observable collisions. We shall focus on cosmic strings, but most of our 
discussion will apply to domain walls as well. In the next Section we specify our assumptions 
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and approximations and find the conditions for a string loop or a spherical domain wall to 
collide with the bubble. In Section III, we study the propagation of the defects that do 
get inside the bubble and find their asymptotic shape. The spatial distribution of defects is 
discussed in Section IV, and the number of potentially observable collision events is estimated 
in Section V. The distribution of collision regions on the observer’s sky is discussed in Section 
VI, although the investigation of the nature of observable effects due to collisions is left for 
future study. Our conclusions are summarized in Section VII. Some technical details of the 
calculations are given in the Appendices. 

II. BUBBLE-DEFECT COLLISIONS 

A. Bubble spacetime 

In this section, we shall discuss the conditions for a string loop (or domain wall) nucleated 
outside of our bubble to collide with the bubble. To set up the problem, we assume the 
spacetime outside the bubble, which we shall refer to as the false vacuum, to be a de Sitter 
space with a constant expansion rate Hp. It can be partially covered by fiat chart metric 

ds^ = dt^ — [dr^ + > (1) 

where dfl^ = dd^ + sin^ Odcf)^. It will also be convenient to use Cartesian coordinates (x,?/, z), 
related to (r, d, 0) by 

dr^ + r^dU^ = dx'^ + dy^ + dz^. (2) 

The spacetime inside the bubble has the geometry of an open FLRW universe, 

ds^ = dr^ — (r) {d^'^ + sinh^ idVL^ . (3) 

We shall refer to it as true vacuum and approximate it by a de Sitter space with an expansion 
rate Hp < Hp, 

a (r) = sinh(idrr). (4) 

The parameter Hp has the meaning of the slow-roll inflation rate inside the bubble at the 
early stages of its evolution. At later times, the form of a(r) should be different, correspond¬ 
ing to radiation and matter eras. However, we shall see that the shape of strings gets frozen 
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in comoving coordinates soon after they enter the bubble, so the late-time behavior of a(r) 
is unimportant for our consideration. 

We shall disregard the gravitational effects of both the bubble wall and the defect. Fur¬ 
thermore, we shall assume that the bubble wall thickness and its initial radius at nucleation 
are both small compared to the de Sitter horizon Hp^. Then the worldsheet of the wall is 
well approximated by the forward light cone, 

r = 1 - (5) 

where we have chosen the nucleation center to be at the the origin of the flat slicing coor¬ 
dinates given by Eq. Q, namely, r = t = 0. 

B. Defect nucleation 

The nucleation rate of strings and walls is 

T oc e~^^, (6) 

where Fe is the Euclidean action of the corresponding instanton. In the case of strings, the 
Euclidean worldsheet in the thin wall limit has the form of a maximal 2-sphere of radius 
Hp^ embedded in a 4-sphere of the same radius, and the instanton action is 

Se = AniiHp'^, (7) 

where /i is the string tension. The string worldsheet after nucleation is obtained as an 
analytic continuation of this instanton^. For a loop centered at a; = on the x-axis, it is 
given by |T2] 

{x — Tnf + y'^ + -I- z = {x — r„) tan a, (8) 

where is a constant parameter related to the nucleation time and a is the angle between 
the plane of the loop and the a:?/-plane. The point with coordinates (t, x, y, z) = (tn, Tn, 0, 0) 
is called the nucleation center of the loop. It has the important property that its future light 
cone asymptotically approaches the worldsheet of the loop at t —)■ oo. 

^ One can show that these are solutions of the Nambu-Goto equations of motion of strings in a de Sitter 
background. 
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The physical radius of the loop at time t is 


R(t) = ^ 


1/2 


(9) 


It asymptotes to Hp^ at t —)■ —oo and expands as R{t) oc exp{Hpt) at t 3> The part of 
the worldsheet at t —)■ —oo is unphysical and should be cut off at some t ~ However, 
unlike the case of bubble nucleation, the loop radius in ([^ changes monotonically with time, 
so there is no natural choice for the moment of nucleation. The proposal of Ref. na is that 
the loop is formed at the time tj such that the classical action, evaluated from t —)■ —oo to 
t/, is S' 1. This gives 

tn-tf Hp^XnSE, ( 10 ) 

where S'^; is the Euclidean action Q. The semiclassical description of string nucleation 
requires that S'^; > 1. On the other hand, we shall be interested in the situation when Se 


is not too large, so that the nucleation rate is not too strongly suppressed. Then Eq. (10) 
suggests that is a few Hubble times earlier than tn- 

Another issue related to the worldsheet cutoff at some time is that such a cutoff selects 
a reference frame where the loop nucleates. Analysis in Ref. [13] suggests that this frame 
is determined, at least in part, by the initial conditions at the onset of inflation.^ The 
existence of a preferred frame breaks the de Sitter invariance of the false vacuum. This 
memory of the initial state extends indefinitely into the future; it is a manifestation of the 
so-called persistence of memory effect [1]. We shall assume that the coordinate system in 
([^ is chosen so that the loops nucleate on constant-t surfaces. 

The number of string loops of radius R, whose centers are located in a 3-volume element 
dV on any constant-time surface, is given by [ 12 ] 

-7 D 

( 11 ) 


dN = 


where we have dehned A = THp^ and we assumed a lower cutoff at R Hp . The shape of 
the distribution near the cutoff depends on the details of the nucleation process, so strictly 


speaking, Eq. (11) applies only for R^H 


-1 


^ The frame of loop nucleation may also be partly determined by the rest frame of the detector used to 
observe the loop [HIIS]. In our case, however, the role of the detector is played by the bubble, which does 
not have a rest frame. The loop collides with different parts of the bubble, moving at different speeds. 
The results of Refs. [ung cannot therefore be directly applied to this case. This issue requires further 
study. In the meantime, here we adopt the simple model outlined in the text. 
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The description of domain wall nucleation is very similar. The Euclidean worldsheet in 
this case is a 3-sphere of radius Hp^ and the instanton action is 


Se = 27rVif 


-3 
F 5 


( 12 ) 


where a is the wall tension. The wall radius as a function of time and the time of nucleation 
tf are given by the same Eqs. 0.0. 

True vacuum bubbles may interfere with defect nucleation. As we mentioned in the Intro¬ 
duction, defect solutions may not even exist inside the bubbles. Here we shall assume that 
defects cannot be formed if they overlap with any bubbles at the time of their formation.^ 
The defect radius at the time of formation is Rf ~ Hp^, and the nucleation process can be 
pictured as quantum tunneling from R = 0 to R = Rf [ 12 ]. We shall assume that defects 
can be formed only if they do not encounter any bubbles in the course of the tunneling - 
that is, if there are no bubbles in the defect interior at t = tf. On the other hand, once a 
defect is formed, true vacuum bubbles can nucleate in its interior. 


C. Conditions for bubble-defect collision 

Nucleating loops of string will have random orientations, with the angle a in ([^ ranging 
from 0 to 7r/2. For simplicity, here we only consider the case of a = 0. In the case of domain 
walls, the wall radius satisfies the same Eq. ([^ as the radius of a string loop. The great 
circle on the wall, which lies in the plane passing through the center of the bubble, evolves 
in exactly the same way as a string loop with an inclination angle a = 0. For dehniteness, 
in this and the following section we shall focus on strings, but the results will also apply to 
domain walls and will actually be more complete for walls (since we only consider the case 
of a = 0). 

In spherical coordinates (t, r, 0,0), the worldsheet of a string loop with a = 0 can be 
written as 

9 2sinY„cos0 cosr/n —cosY„ ^ ^ , / n 

^2-An-^- in -^ ^ 9 = 7r/2, 13 

cos rjn + cos Xn COS Tjn + COS Xn 

^ Even if string solutions do not exist in true vacuum, open string segments could still nucleate outside the 
bubble with their ends attached to the bubble wall. Similarly, strings that can exist inside but not outside 
the bubble can nucleate with their ends attached to the bubble wall on the interior side. Nucleation and 
evolution of such defects will be discussed elsewhere. 
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where, for later convenience, we have re-expressed the coordinates tn and of the nucleation 
center in terms of its conformal coordinates in the closed slicing of de Sitter space dehned 
by 


^ — Hptn _ 


sm rjn 


r„ = 


cos rin + cos Xn 
sin Xn 


^ ( 14 ) 

cos rjn + cos Xn 

with 0 < < TT and —tt < < 0. Assuming that the nucleation center is in the part of 

de Sitter space covered by the flat-slicing coordinates Q, we must have Xn ^ Vn + tt. All 
coordinate systems that we use in this paper and the relations between them are summarized 
in Appendix A. 

To derive the conditions for a string loop to cross the bubble wall ([^, we shall focus on 
the x-axis, y = z = 0, where we can And the nearest point and the most distant point on 
the loop with respect to the bubble center. The intersections of this axis with the string 
loop and with the bubble are 


sm 


cos r]n + cos Xn 
xb± (t) = ± (l - , 


± 4 / + 


sm r]n 


{cos r]n + cos Xn) 


2 ’ 


( 15 ) 

( 16 ) 


where the subscript S denotes the string loop, B denotes the bubble, and ± denote the right 
and left parts of the string loop/bubble respectively. For definiteness, we shall assume that 
the loop nucleation center lies on the positive side of the a;-axis, that is, at x = > 0. 

Then, in order to have a collision, we need 


I xs- {t —)■ cx)) I < 1, (17) 

or —TT 12 < Xn + Vn < 7r/2, and the string loop will completely enter the bubble if it also 
satisfles"^ 


Xs+ {t ^ oo) < 1, (18) 

or Xn — ^n < 'x/2. The collision starts on the right side of the bubble (that is, on the positive 
half of the x-axis) if 

xs- (0) > 0, (19) 

^ Since the loop nucleation center is assumed to be at a; > 0, it is obvious that we always have xs+ > 0. 


7 









FIG. 1: The different regions of the parameter space of {tn,rn) shown over the conformal diagram 
of de Sitter space. String loops with nucleation centers in the blue and green regions will partially 
enter the bubble, and string loops with nucleation center in the red region will completely enter 
the bubble. String loops with nucleation centers in the green and red regions will hit the bubble 
from the left, and string loops with nucleation center in the blue region will hit it from the right. 


or Xn > 7r/2; otherwise, the collision, if there is one, will start on the left side of the bubble. 
Finally, by solving xs- = xb±, we hnd the time of the collision. 


tc± = 7 ^ log 
Hf 


1 + 


cos Xn 


( 20 ) 


cos77„ =Fsinx„_ 

There is always a solution of xs- = xb±- However, for tc± to have physical meaning, we 
need tc+ > 0 if the collision is from the left or tc- > 0 if the collision is from the right, since 


the bubble worldsheet described by Eq. (16) with t < 0 is unphysical. 

To summarize, we project all the above constraints into the parameter space of (//„, Xn) 
in Fig. The bubble nucleation center is at {r]n,Xn) = (—7r/2,0), and the bubble wall is 
approximated by the future light cone of this point. String loops whose nucleation centers 
are inside the bubble (that is, in the red area in the diagram) will completely enter the 
bubble. This follows from the fact that the future light cones of such nucleation centers 
are completely contained within the bubble. On the other hand, the future light cones of 







nucleation centers in the green and bine areas are only partially contained in the bnbble. 
The corresponding loops will therefore only partially enter the bubble. String loops with 
nucleation centers in the blue area of the diagram hit the bubble from the right, while loops 
with nucleation centers in the red and green areas hit it from the left. In the latter case, 
the loop must encircle the bubble; hence, it has to be formed at f < 0, prior to bubble 
nucleation. 

We emphasize that loops can form in the false vacuum even when their nucleation centers 
are in the true vacuum inside the bubble. For example, the loop nucleation center can be 
in the red area of the diagram, while the loop formation time, which is earlier than can 
be at < 0 (that is, rjf < —7i/2). The loop is then formed prior to the bubble nucleation. 
When the bubble forms, it is initially encircled by the loop, and as it expands, the loop is 
eventually completely engulfed by the bubble. 

We note that the division of spacetime into different regions in Fig. is invariant under a 
class of de Sitter boosts. This is obvious for the boundary between the red and green regions, 
which corresponds to the bubble wall, but may need a demonstration for the boundary 


separating the blue and green regions. According to Eq. (19), this boundary is given by 
Xn = or 

Tn = -n/I + _ (^21) 

Using the embedding coordinates (X, U, Z, hF, U) dehned in Appendix A, this can be rewrit¬ 
ten as 


W = 0. (22) 

On the other hand, the bubble wall lies at IF = Hp^. de Sitter transformations that leave 
the wall invariant correspond to bD Lorentz boosts in the X, Y and Z directions, 

^ - PX, 

IF' = IF, (23) 

where (3 is the boost constant. These transformations do not affect the IF coordinate and 
thus leave the boundary ( [2^ invariant. (Boosts in the IF direction change the nucleation 
time parameter of the loop.) 
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The boost invariance of the boundary between blue and green regions can also be under¬ 
stood geometrically. Defects with nucleation centers in blue and green regions hit the bubble 
respectively from the right and from the left. By continuity, defects with nucleation centers 
on the boundary between these regions must hit the bubble at the moment of its formation. 
Worldsheets of such defects should pass through the origin of the light cone representing the 
bubble. This is a geometric fact which is invariant under boosts. 


III. STRING PROPAGATION INSIDE THE BUBBLE 

We would now like to study the dynamics of the string after it enters the bubble. We 
assume that the fields making up the string do not interact with the field of the bubble, so 
the string responds only to the change of spacetime geometry. We shall start with the case 
of Hx = Hp, and then extend the analysis to the case of Hp < Hp. 


A. Hp = Hp 


In this case, the whole spacetime is a pure de Sitter space; hence the worldsheet of the 


string loop can be described by Eq. (13) both inside and outside the bubble. The only 


difference is that the observer inside the bubble uses the open-slicing coordinates, which 
relate to the flat coordinates used earlier by 

1 


cosliTTr = coshift- 

2 


sinh H T sinh ^ = e 




(24) 


where Hp = Hp = H. Thus, according to the observer inside the bubble, the worldsheet of 


the string found earlier in Eq. (13) takes the form 


cos 7]n cosh ^ tanh Hr — sin Xn cos (j) sinh ^ tanh Hr — cos Xn = 0. 


(25) 


As illustrated in Fig. the shape of the string can be found by setting r = const in this 
equation. At r S> H~^, tanh Hr ^ 1, and thus the string freezes in the comoving coordinates 
within a few Hubble times after it enters the bubble, approaching the asymptotic shape 


cos rjn cosh ^ — sin Xn cos 0 sinh ^ — cos Xn = 0- 


(26) 
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FIG. 2: Intersecting worldsheets of the bubble wall and of a nucleated string loop shown in the 
conformal coordinates. The pink surface is the worldsheet of the bubble, the green surface is a 
constant time surface in the open-slicing coordinate system inside the bubble, and the blue surface 
is the worldsheet of the string loop. The size of the string loop in conformal coordinates shrinks with 
time, but the actual physical size grows with time. Highlighted by the red line is the intersection 
between the green and blue surfaces, showing the shape of the string at a certain time. The string 
has infinite length and a fixed asymptotic angle at spatial infinity. 


Eq.(26) has two solutions for 


g«± = 


2 cos Xn \/l — 2 cos 2r]n + cos 2(j) + 2 cos 2xn sin^ (j) 
2 cos rjn — 2 cos (j) sin Xn ’ 


(27) 


which are related by (0) = —(0 -1- tt); the physical solution is chosen by requiring ^ > 0. 
Thus, strings nucleating in the blue region are described by and strings nucleating in 
the green region are described by For some of the strings nucleating in the red region, 
we need the positive parts of both and 

We can define the minimal distance from the string to the origin as at 0 = 0: 




log 


'cosXn + smr]n 


e^tn -f- r„) — 1 

cos rjn - smXn_ 

log 

eHt„ (^1 _ 


(28) 


where we also give the expression of in terms of and for later use. According to the 
relation between and a negative indicates that the minimal distance appears on 
the left half of the x-axis. To clarify the meaning of it is helpful to project the strings 
at r —)• cx) onto the Poincare disk. Introducing zp = tanh | and (j)p = (j), we can project the 
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FIG. 3: The shape of strings at r —oo projected on the Poincare disk. We plotted the shape of 
3 different kinds of loops and labeled the corresponding and Acj). The parameters {'r]n,Xn) are 
(—7r/4,7r/2) for the blue loop, (—7r/8,7r/4) for the red loop and (—Svr/S, 7r/4) for the green loop. 


plane of 6^ = 7r/2 at constant time r onto a nnit disk. Then a string at r —)■ cxd has the form 


zp = zq cos (pp + \JRp — Zq sin^ 0p, 


(29) 


which is a circle on the Poincare disk. Here, we have dehned zq = 


smxn 


Rp = 


sin ri„ 




. In Fig. 


, = and 

COSXn+COS?7n “ 

we show some examples of different kinds of loops 


cos Xn+cos ri„ 

and mark the corresponding distances The asymptotic shape of the loop depends on 
the location of its nncleation center {jj^^Xn)- 

For a loop nncleated in green or bine regions in Fig. the string will appear to have 
inhnite length, with its two asymptotes becoming straight and spanning a hxed angle from 
the origin. This angle can be fonnd by sending ^ cxd, 

’(1 + rD - 


A(f) = 2 arccos 


cos?7„ 
sin Xn 


= 2 arccos 


2rr. 


(30) 


For a nncleation center in the red region of Fig. [T| the loop initially snrronnds the origin, 
then shrinks and hnally freezes in comoving coordinates, bnt physically always expands. At 
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FIG. 4: An example of string worldsheet with nucleation center in the blue region of the conformal 
diagram (Fig. in the case when Ht = Hp. The thick blue and black solid lines are the worlsheets 
of the string loop and of the bubble wall, respectively. The loop is formed outside of the bubble and 
then partially enters the bubble. The blue dashed lines denote the light cone of the loop nucleation 
center, which is represented by a blue dot. The black dashed lines indicate the surfaces of t = 0 
and t = —oo of the flat slicing. Finally, the dot-dashed lines show the boundaries between different 
parameter space regions of Fig. 


r —)■ oo, we can find the maximal distance from the origin to the loop: 




-log 


'cosy^ -bsinr;^' 


(1 — Tn) — 1 

cost^n + siny„_ 

log 



( 31 ) 


Some examples of loop worldsheets with different nncleation centers are shown in Figs. 


B. Ht < Hp 


In this case, we cannot nse Eq. (13) after the string enters the bubble. In order to extend 


the string worldsheet into the bubble, we use our assumptions that the bubble wall gravity 
can be neglected and that the initial size of the bubble is small compared to the horizon. 
Furthermore, the bubble interior is initially curvature dominated, so a(r) ~ r, independent 


® We use a 1 -|- 1 dimensional version of the conformal diagram for de Sitter in these figures in order to be 
able to illustrate both sides of the string worldsheet. 
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FIG. 5: An example of string worldsheet with nucleation center in the green region of the conformal 
diagram in the case when Ht = Hp. The loop initially encircles the bubble and then partially 
enters the bubble. 


of H. We can therefore extend the false vacnum spacetime ontside the bubble into the bubble 
interior up to some small time The string shape and velocity at 0 < Tj < r^, 


namely ^ and ^ (Tj,0), can then be found from Eq.(13). Hence, as —)• 0, 

should satisfy 


A(l]nyXnytlp)riCOSh^-B(T]nyXnyHp)TiCOSii)Smh^, ^ 1 , 


(32) 


where we have defined 

A{r]n,Xn,H) = H -, 

cos Xn 

B {rin,Xn,H) = HtaiaXn- (33) 


With these initial conditions, the following evolution inside the bubble can be found by 
solving numerically the equation of motion of the string. 

The solution can also be found analytically as follows. We already know that for a string 
moving in a pure de Sitter space with H = Hp, there is a class of solutions which satisfy 


cos cosh ^ tanh Hpr — sin Xn cos 0 sinh ^ tanh Hpr — cos Xn = 0- 


(34) 


Eq. (34) is the same as Eq. (25), except for a difference in H. We then observe that the 


solution described by Eq. (34) will have the same shape and velocity as the string described 
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FIG. 6: An example of string worldsheet with nucleation center in the red region of the conformal 
diagram in the case when Ht = Hp. The loop initially encircles the bubble and ends up completely 
inside the bubble. 


by Eq. (32) in the limit of r —)■ 0, as long as 


A {T]n, Xm Hf) = A Ht) , 

B {r]n, Xn, Hf) = B (77^, x'„, Ht) . 


(35) 


In other words, for a loop nucleated air] = rjn and y = x„ in the false vacuum, the evolution 


of its part inside the bubble can be described by Eq. (34), with t]'^ and x'n found through 


Eq. (35). This gives 


COSr]n COS77„ , 

= 7-r tan Xn = 7 tan x„, 


cos Xn cos xh 

where 7 = Ht/Hf- We also express this mapping in terms of and for later use: 

, 2rnC^ 


c' = 


(1 - 7)r2c2 + (1 + 7 )c2 + 7 - 1 ’ 

_ (1 -7)cV; +c^(l +7) -1 + 7 _ 

^ (7'^ — 1) (r^ — 1)^ + 2 (1 — 7^) + 2 (1 + 7'^) — (1 — 7^) 


(36) 


(37) 


where we have defined c = and c' = For 7 = 1, we get back to the case of 

Vn = Vn and Xn = Xn- For 7 -)■ 0, Xn approaches 7r/2. 

We have verified that this mapping between solutions agrees well with a numerical cal¬ 
culation. We note, however, that we cannot find a corresponding (?7(j,Xn) for any {rin,Xn)- 
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Xn Xn 


FIG. 7: The parameter space of {i]n,Xn) with 7 = 0.5, 0.1. The color coding is the same as in 
Fig.[Tl except that we added a gray color, marking the nucleation centers of loops which collapse 
after they enter the bubble. 



FIG. 8: A string worldsheet with the same nucleation center as that shown in Fig. [^in the case 
when 7 = 1/2. 
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FIG. 9: A string worldsheet with the same nucleation center as that in Fig. in the case when 
7 = 1 / 2 . 



FIG. 10; A string worldsheet with the same nucleation center as that in Fig. in the case when 
7 = 1/2. We do not show the nucleation center in this case, because it is in the fictitious false- 
vacuum de Sitter space. 
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For example, for Xn = 0 we have Xn = 0; Eq. (36) gives cos rjn = 'y cos rj'^. This has 


no solutions for r]'^ when cos rjn > 7 . More generally, it can be seen from Eq. (35) that, for 
iVn, Xn) such that 


cos + tan^ Xn cos Xn, 


( 38 ) 


there is no corresponding {'f]n,Xn)- Loops with cos r]n = + tan^ Xn cos Xn have r]n = 0] 

their physical radius shrinks to the horizon size, at r —>■ 00 . And loops satisfying 


(38) collapse to a point in a hnite time. Such loops enter the bubble completely with radii 
smaller than and typically collapse on a timescale < - unless the inequality (38) 


is only marginally satished. The parameter values corresponding to loops collapsing after 
they enter the bubble are marked by gray color in Fig. Comparing the diagrams with 
7 = 0.5 and 7 = 0.1, we see that, as 7 decreases, more and more red regions turn into gray. 
This is because more and more incoming string loops tend to collapse as the horizon inside 
the bubble increases. We have verihed that for the rest of the parameter values {rjn, Xn) that 
lead to a collision, one can always hnd the corresponding {rjn ^ Xn')- 

In Figs. |8 IT, we show the worldsheets of strings with the same nucleation centers as in 
Figs. |4]j^ continued into a bubble with 7 = 1/2. Note that the slope of the string worldsheet 
appears to be discontinuous at the bubble crossing. This is a spurious effect, which is due 
to a discontinuous derivative (in the direction normal to the wall) of the conformal factor 
that we used to generate the conformal diagram. 


We note hnally that the mapping prescription (36) was derived only for the case of 


a = 0. We were not able to generalize it for a nonzero a. This method can therefore always 
be applied to domain walls, but propagation of strings with a 7 ^ 0 can only be studied 
numerically. 


IV. DISTRIBUTION OF DEFECTS 

Having in mind possible observational effects, we shall now discuss the statistical distri¬ 
bution of strings and walls as it appears to an observer inside the bubble. The potentially 
observable effects can be roughly divided into two kinds. If the defect is inside our Hubble 
radius, it can be observed directly, for example through gravitational lensing effect. In this 
case, one may be interested in the spatial distribution of strings or walls at a given time r. 
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Another possibility is that the defect is no longer within our Hubble radius, but we are in 
the future light cone of its collision event. We may then hnd some observational signals of 
the collision, for example, a distortion of the spectrum of density perturbations, or a burst 
of gravitational waves emitted during the collision. Then the number of accessible events is 
just the number of collisions within the past light cone of the observer. 


A. Infinite domain walls 

From the point of view of an interior observer, string loops and domain walls that partially 
enter the bubble can be characterized by the parameters and A(^. In the case of strings, 
there is also a third parameter, a, which is the angle between the plane of the loop and the 
line connecting the loop and bubble nucleation centers. We shall hrst consider the simpler 
case of domain walls. Since the walls get frozen in comoving coordinates shortly after they 
enter the bubble, their distribution at a late time r, e.g., at the end of inflation, is well 
approximated by the distribution at r —)■ cxd. 

We start with the distribution of domain walls with nucleation centers at {tn,rn)- Then, 


using Eqs. (28) and (30) with H = Ht, we can hnd the distribution with respect to ‘^4’)- 


The calculation is somewhat tedious but straightforward. The number of domain walls with 
nucleation centers in spacetime volume element 


dQ,n = e^^‘'^"A'Kr‘^drndtn 


(39) 


IS 


dN = 


( 40 ) 


where T is the nucleation rate per unit spacetime volume. Then we have 


dN = THp^e^^^^^Anrldrndtn 
= A-K\c?r^drndc 




d (r„, c) 


d^md^fpi 


(41) 

(42) 


where, as before we have dehned 


A = 


(43) 
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FIG. 11: The distribution of and A(f> for infinite strings with two different values of 7 . The 

distribution is only defined in the region below the red dashed line. The numbers labeled on the 

, , dN/X 

contours are log 


Using the Jacobian calculated in Appendix we finally obtain 

dN „„ , 7 ^ (2cosh,^mSin ^ — sin A0sinh,^m) 


= 327rA- 


d^mdAcj) (2 cos ^ cosh^m — 2sinh^m)^ + 27^(1 — cos A0) 

This distribution is plotted in Fig. It is dehned only in the region 


5/2 • 


(44) 


7cos^ + U1-(1-72)cos2^ 

- - -at- > 

1 - cos ^ 


(45) 


which is below the dashed red line in the hgure. This is because the fiat chart of de Sitter 
space covers only the region of Xn — Vn ^ tt. 


At a hrst glance, Eq. (44) may suggest that the number of defects is strongly suppressed 
by a factor of 7 ^ for small 7 . However, as can be seen in Fig. 11, there is a band near 


tanh,^m = cos A0/2 where the distribution is enhanced. In the vicinity of this band. 


dN 


2ti\ 


COsh^ 


(46) 


d^rndAcj) 7 

A useful characteristic of the distribution is the total number of defects N{^) in the region 
within a distance ^ from the origin. This is given by 




d^n 


"2- ^ dN 

dA(j)- 


dvrA 


'A 


'0 


d^mdAcj) 


[2e+(l + 7') Sinh2e]. 


(47) 


20 



























Here, the factor 1/2 accounts for the fact that the distribution (44) is dehned only in a half 
of the full parameter space (A0, 

Observational bounds on the curvature of the universe indicate that our present Hubble 
radius cannot exceed 10% of the curvature radius. This means that a central observer can 
detect only defects located a-t ^ ^ 0.1. The number of domain walls; within the 

observable region should then satisfy 




( 48 ) 


With A < 1, we thus expect no more than 0{1) walls in the observable universe. This 
number may be somewhat enhanced if there is a large number of different kinds of domain 
walls, as in the axiverse picture [TB] . 

We note hnally that the analysis in this section does not account for the constraint that 
the defects should be completely outside the bubble at the time of their formation. Some of 


the domain walls included in the distributions (47) and (48) do not satisfy this constraint. 


and thus these distributions overestimate the number of defects. We expect, however, that 
they still give the right order of magnitude estimates. (We have verihed that this is indeed 
the case for 7 -C 1.) 

B. Infinite strings 


The calculation of domain wall distribution in the preceding subsection cannot be ex¬ 
tended to the case of strings because, as we mentioned at the end of Sec. HI.B, the mapping 


prescription (36) can only be used for strings with inclination angle a = 0. In order to see 


the qualitative differences introduced by a variable a, we shall restrict the analysis to the 
case of 7 = 1, discussed in Sec. HI.A, where such a mapping is unnecessary. 

The distribution of string loops in the parent vacuum is given by 


dN = dvrAr^m dc (i(sina). 


( 49 ) 


where c = as before. A loop with a nucleation center at {tn,Tn) and inclination angle 

a will overlap with a sphere of radius r centered at the origin if a < am{rn, c, r), where 


am = arccos 




( 50 ) 
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The total number of loops within a comoving radius r from the origin is thus given by 

N{r) =ATrX / dcc^ sinam(?"n, c, r), (51) 

J D 

where the integration domain D is specihed by the conditions 

(52) 

(53) 


c ^ - r <rn<c ^ + r, 


> 1 - c \ 


c > 0 . 


(54) 


The hrst of these conditions requires that a loop with a nucleation center at (t„, r„) and 
a = 0 overlaps with the sphere at t —)■ oo. Some loops with nonzero a will not overlap with 
the sphere; this is accounted for by the factor sinctm- The second condition requires that 
the loop nucleation center is outside of the bubble. 

The integrations in Eq. (51) are rather tedious. The details are given in Appendix [C| 
here we only give the hnal answer: 

A 7 r^r^(r^ + 16r + 4) 


N{r) = 


(56) 


In the asymptotic region t —)■ cx), the comoving radius r can be expressed in terms of the 
^-coordinate inside the bubble, 


r = 


sinh^ 

1 + cosh ^ 


Substituting this in (55), we obtain 

Avr^ 


sinh^ ^ (3 -|- 5 cosh 16 sinh . 


(56) 


(57) 


An important difference of this distribution from Eq. (48) for domain walls is its behavior 
at small 

iV(0 ^ (e « !)• (58) 


The expected number of strings within the observable region is thus 


iV(a) < O.IA. (59) 

This is small for any A < 1. 

The difference between the wall and string distributions at small ^ is due to the fact that 
string loops with a sufficiently large inclination angle a do not cross a small sphere near the 
origin, even if a domain wall of the same radius would cross it. We expect that the string 
distribution for 7 < 1 will exhibit the same qualitative behavior. 
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C. Closed strings and walls 


Closed strings and walls, which completely enter the bubble, have their nucleation centers 
in the bubble interior. We begin by considering a loop or domain wall whose nucleation 
center is at = 0. The corresponding worldsheet is given by 

1 


cosh^ = 


tanh Htt tanh Htt' ’ 


where 


( 60 ) 


tanh Htt^ = 7 ^ tanh HpTn. 


The size of the defect can be characterized by the radius 


( 61 ) 


R = Hrp ^ sinh Htt sinh ^ = Hj, ^ cosh Htt (coth^ HtTu — coth^ Htt) . 


(62) 


The worldsheet of a defect with a nucleation center at {Tmin) can be obtained from (60) 


by a de Sitter boost (23) with (3 = tanh^„. The resulting loop or domain wall will look. 


respectively, like a circle or a sphere of radius (62) from the new nucleation point 


As in the preceding subsection, we shall disregard the constraint that defects can only 
form outside of the bubble. Our result will therefore be an overestimate of the density of 
closed defects. We shall see, however, that this estimate still tends to be rather small. 

The number of defects with nucleation centers in the interval (dr^, d^) is given by 


dN = rdf24 = ATiXHpsirhr ^d^sirhi HpTndTn 


( 63 ) 


Note that the 4-volume element dVL^^ is in a hctitious de Sitter space, obtained by continu¬ 
ing the parent false vacuum into the bubble interior. The nucleation centers are uniformly 
distributed in this extended region (with a constant density per unit 4-volume). This dis¬ 


tribution is invariant under de Sitter boosts (23), which implies that the defect distribution 


should have the same symmetry. With these assumptions, the defects are uniformly dis¬ 
tributed on the surfaces of constant time r inside the bubble.® 


® The cutoff at t/ would amount to requiring that defects can only be formed with t/ < 0, or equivalently 
with tn < Se- This cutoff would of course break the de Sitter invariance. We expect that it would 
suppress the closed defect density at large values of The effect of the cutoff at tj on the number of 
observable collisions will be discussed in Section V. 
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To find the distribution of defects of a given radius i? on a hypersurface of constant r, 
we express r^, in terms of R from Eq. (62), 


dN = 


Xj^dV 


RnRdR 


( 64 ) 


tanh* HtTs/W+M ’ 

Here, dV = Ansinh^ Htt s'mh^ ^d^ is the volume element on the hypersurface, and 
Ro{t) = smhHTT is the curvature radius of the spacelike slice r = const. For 7 1 

and Htt 3> 1 , the distribution simplihes to 

RoRdR 


dN ~ X-f^dV- 


(65) 


We can also dehne the comoving size as R = R/Ro and the comoving volume dV = 
dV sinh“^ Htt. Then 

RdR 


dN ~ X-f^dV- 


( 66 ) 


{R^ + 1)3 

which simply says that the defects are formed at early times r and are then stretched and 
diluted by the expansion, leaving the comoving distribution unchanged. 


The distribution (64) shows that the number of defects is strongly suppressed by 7 *^. This 
suppression is in qualitative agreement with the shrinking of the red region in Fig. Hence, 
we should not expect to see any closed loops or walls from the multiverse within our horizon, 
unless 7 rsj ^ 

Closed strings and walls can also nucleate in the true vacuum inside the bubble. The 
corresponding nucleation rates are given by Eq. ([^ with the instanton action ([^ or (12) and 
with Hp replaced by Hp. This rate, however, will be exponentially suppressed for Hp < Hp, 
assuming that the string or wall tension is the same inside and outside the bubble. If the 
true vacuum admits lighter defects, which do not exist in the false vacuum, they could have 
a higher nucleation rate and could produce some observable effects (see e.g. CZI)- However, 
observation of such defect would not provide evidence for a multiverse, so we do not consider 
them here. 


V. THE NUMBER OF COLLISION EVENTS 

We shall now estimate the number of collision events within the observer’s past light 
cone. We consider a central observer, located at .^ = 0 inside the bubble. The past light 
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FIG. 12: The red cone is the past light cone of the observer which intersects the bubble wall on 
the sphere S, and the blue plane is the hypersurface S. Black circles show the domain walls on 
that hypersurface. 


cone of such an observer, originating at time tq, is given by 


g _ tanh(ifr'7'o/2) 
tanh(i?rr/2) 


( 67 ) 


Its physical radius at time r < tq is 

R{t) = sinh^ sinh H^t = coth.{H'rTo/2) 

The bubble wall in these coordinates is at r = 0; it intersects the past light cone on a sphere 
S of radius 


cosh {Htt/2) 
cosh^(HTTo/2) 


( 68 ) 


R{0) = ta.iah{HTTo/2). (69) 

For To 3> we have i?(0) ~ - which is to be expected, since is the horizon 

radius in the de Sitter space of the bubble. 

Let us consider the particular hypersurface of constant time S : f = const in the flat 
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slicing coordinates given by Eq. Q, which includes the sphere S. From Eq. ([^ we have 

1 + 7 


e^Fti ^ 


7 


(70) 


This hypersurface extends from the sphere into the false vacuum (see Fig. 12). Now, the 


number of defects of radius R, whose centers are located in a 3-volume element dV on any 
constant-time surface, is given by Eq. 01. 


dN = A^dV, 


(71) 


with a lower cutoff at i? ~ Hp^. The defects that intersect the sphere S at t = ti must have 
collided with the bubble wall at some earlier time tc- We shall estimate the number of such 
collisions, focussing Erst on domain walls. 

It will be convenient to use the units in which Hp = 1 and Hp = 7. The collision time 
tc for a domain wall with nucleation center at {tn,rn) can then be found from 

2^ 


^tc tn 


l-z^' 


where z = (r„ — l)e*". It should satisfy 


tf ^ tc ^ tp 


(72) 


(73) 


where the time of wall formation t/ is related to tn by Eq. (|10[). With the aid of Eqs. (72), 


(10), the condition tc> tf can be expressed as 


ic if _ in if 


2z 


1 — z‘ 


;Se > 1 , 


or 


For S'e S> 1, this gives 


^ \ S% + 1 - 


l)e'" > 


2Sf 


(74) 


(75) 


(76) 


The parameters tn, can be expressed in terms of the domain wall radius R and the physical 
distance Vphy from its center to the bubble center,^ 


e *" = 


1 + 7 


vr^ - 1 , 


(77) 


^ Note that the distance rphy is calculated in the fictitious de Sitter space, obtained by extending the false 
vacuum to the bubble interior. 
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6 ^phy 


7 


1 + 7 


^phy 


Substituting this in (76) we obtain 


'l^phy > h 1 + —— 1. 

7 


(78) 


(79) 


For the most numerous domain walls with radii near the lower cutoff, R Se and the 
above condition simplihes to 


^phy ^ ^ S” 

Similarly, requiring that < ti, we obtain the condition 


< - + R. 
7 


(80) 


(81) 


It simply states that the domain wall should be close enough to overlap with the bubble. 

The total number of intersections can now be found by counting all domain walls in the 
distribution (0 that satisfy the conditions ([^, ([^, 


A7> A 


dR 


n-'^+R 


27r 


— / _ dvrrphydrphy ~ — Ay 


-2 


(82) 


'7 ^ + 1 


where in the last step we assumed that 7 1. We see that the number of potentially 

observable collisions is enhanced by a factor of 7 “^ compared to the number of walls (48) 
that are physically present within the visible universe. A similar enhancement factor was 
found in Ref. [7] for observable bubble collisions. For 7 +; 1 it is possible to have A/" S> 1, 
provided that the wall nucleation rate A is not very strongly suppressed. 

The total number of intersections for cosmic strings TVs can be found by a calculation 


similar to that in Sec. IV.B. One obtains an integral similar to Eq. (82), but with a weighting 
factor sin0^(7?, ''"phy, 7 ); due to the variable inclination angle, am is the maximal inclination 
angle for given R and rphy. 


r.v,2 


Urn = arccos 


phy 


+ R2-7 


- 2 n 


2^ phy 7? 


(83) 


For small 7 , we hnd 


A/; 


~T^ 


Ay 


-2 


(84) 


In the above analysis we disregarded defects that have their nucleation centers inside the 
bubble. Such defects encircle the bubble at t = 0, and we shall assume as before that they 
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can only be formed at tf < 0. Then their nncleation centers must be at tn < ^ InS'e. 

The number of such defects can be estimated as 


Af' ~ rf] 4 , 


(85) 


where is the 4-volume of the relevant spacetime region, 


ff4 


dr ■ dvrr^ 



Here, 


hir) = — ln(l — r) 


is the worldsheet of the bubble wall and 


( 86 ) 


(87) 


r. ~ 1 - ^ (88) 

OE 

is the comoving radius of the bubble wall at time f*. Performing the integrations and 
assuming that Se 3> 1 , we obtain 

V' ~ yASi. (89) 

The factor S'!; 3> 1 appears to enhance the number of collisions; however, the nncleation 
rate A ~ is suppressed at large Se, so we always have Af' < 1. 


VI. COLLISION MARKS IN THE SKY 


We shall now try to see what the collision events might look like in the sky. We shall 
not attempt any realistic analysis of observational effects that the collisions may produce in 
the CMB. Instead, we shall assume that collisions with defects leave some sort of marks, or 
‘scars’ at very early times in the bubble (r —)■ 0). One can assume, for example, that defects 
disintegrate upon colliding with the bubble and their energy is deposited in close vicinity 
of the impact. Assuming that this energy excess somehow becomes observable, we shall be 
interested in the shape of the resulting ‘scars’ and in their distribution in the observer’s sky. 


The mark left by a collision with a domain wall has the form of a disc. In Fig. 13 we 


show the distribution of wall collision marks in the sky for two values of the parameter 7 , 
using the Mollweide projection of the celestial sphere. Note that this projection distorts 
the circular shape of the marks. The distributions were generated for a hxed number of 



y=0.1 Se=3 N=20 


7=0.01 Se=3 N=100 




FIG. 13: Collision marks of domain walls on the sky. The left panel shows 40 collision events with 
7 = 0.1, and the right panel shows 100 collisions with 7 = 0.01. The corresponding nucleation 
rates are A = 0.025 and A = 0.0016, respectively. 

y=0.1 Se=3 N=20 7=0.01 Se=3 N=100 




FIG. 14: Collision marks of cosmic strings on the sky. The numbers of collision events, the values 


of 7 and the color code are the same as in Fig. 13, Note that all strings in these simulations 
nucleate in the blue region of the parameter space. 

collisions A/”, with a random choice of nncleation centers. We only inclnded walls that are 
completely ontside the bnbble at the time of their formation tf. Bine and green cnrves in 
the fignre are respectively the bonndaries of collision marks dne to walls with nncleation 
centers in the bine and green areas of the parameter space of Fig. The green collision 
marks correspond to domain walls that formed prior to the bnbble formation and initially 
enclosed the bnbble. A typical bnbble can be expected to collide with snch walls only if the 


wall nncleation rate is relatively high, as in the left panel of Fig. 13 Bine collision marks 
are dne to walls that formed ontside the bnbble. 

For bine cnrves the area of impact on the sky is enclosed by the cnrve; in other words, 
the affected area is always less than half of the sky. For green cnrves the affected area can 
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either be inside or outside, but we find that it is predominantly outside. For example, all of 


the green curves in the left panel of Fig. 13 have the affected region outside. 

For a head-on {a = 0) collision with a string loop, the mark is just a straight segment, 
while for a collision with a ^ 0 the mark segment is curved. Two simulated distributions 
of string collision marks in the sky are shown in Fig. As for the walls, the distributions 
were generated for a fixed number of collisions, but now the nucleation centers and loop 
inclination angles a had to be chosen at random. The details of the simulations are given 
in Appendix D. 

The distribution of angular sizes of domain wall collision marks on the sky can be found 
analytically by calculating the 4-volume of the appropriate nucleation region. The result is 
shown in Fig. [^by the red solid line. The calculation involves some tedious integrations, 
and we do not reproduce it here. Instead, we shall find the form of the distribution for 
relatively large angles, which can be calculated as follows. 

As in the preceding section, it will be convenient to use units in which Hp = 1 and 
Hp = 7 . Once again, we consider the spherical surface 5, defined by the intersection of the 
observer’s past light cone and the bubble worldsheet, and the fiat hypersurface T, : t = ti 
that includes S. Now consider a domain wall of radius R ^ 7 “^ lying in the hypersurface 
S with its center at a distance p < R from S. The intersection of the wall with the surface 
S will be seen from the origin at an angle 


A(j) 27 . 


(90) 


Using the wall distribution (0 and expressing p in terms of A0 from Eq. ( [90| ), we have 

dR AcpdAcp 


dN = 27rA7 


-4 




(91) 


To find the distribution of angular sizes A0, we need to integrate Eq. (91) over i?, with 
a lower limit of integration at Rmin = A 0 / 27 . This gives 

327rA dA(j) 


dN = 


(92) 


3 

This should be accurate, as long as 7 “^ 3> Rmin S> 1, that is, 1 S> A0 3> 27 . 

Numerically generated angular size distributions for collision marks made by domain walls 


and cosmic strings are also shown in Fig. lA They were obtained by averaging over 300,000 
instances of randomly chosen defect nucleation centers and using the parameter values of 
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FIG. 15: Numerically generated angular size distributions for 300,000 collision marks made by 
domain walls (the left plot) and by cosmic strings (the right plot). In the left plot, the analytically 
calculated distribution for domain walls is shown by a solid red line. We fit the distributions with 


a power law (93) in the range of 0.01 < Afp < 0.05. (The points used in the fits are shown in red.) 
The best fits (shown by dashed black lines) yield the values p = —3.03192 for domain walls and 
p = —3.00792 for strings. 


7 = 0.001 and Se = 3. The vertical axis in the hgure shows the quantity dN/d{A(j)) (up to 
an arbitrary factor). For domain walls, the numerically generated distribution agrees very 
well with the analytic result. For cosmic strings, we dehned the angular size Acp as the angle 
between the end points of the string collision mark on the sky. In this case, we did not 
attempt to derive the distribution analytically, so we show only the numerically generated 
distribution. 

We htted the distributions with a power law 


dN oc {A^yd{AcP) 


(93) 


in the range® of 0.01 < A(p < 0.05. The best £t yields the values p = —3.03192 for 
domain walls and p = —3.00792 for strings, suggesting that the string distribution is likely 


to approach the same asymptotic form (92). 

The analysis of bubble collisions in Ref. has shown that the distribution of collision 


To reduce noise in the numerical data, we used logarithmic bins for Ap. The fitting range for the power 
law was chosen so that the number of collision marks per bin is no less than 18 for all bins in the range. 
The lower bound of the range, Acp > 0.01, was chosen to satisfy the condition of validity of the power law 
approximation, Acp ^ = 0.002. 
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events on observer’s sky is generally anisotropic. Assuming that inflation begins on a surface 
t = ti in the flat slicing chart ([^, only the central observer, who is comoving with respect to 
the geodesic congruence orthogonal to this initial surface will see an isotropic distribution. 
This effect, which persists even in the limit —>■ cx), has been dubbed “the persistence of 
memory”. If the coordinates inside the bubble are chosen so that the central observer is at 
the origin, ^ = 0, then for all other observers the distribution will be peaked in the direction 
away from the origin. We can, of course, perform a de Sitter boost that brings any given 
observer to ^ = 0, but this transformation will also affect the initial surface, which will now 
slope down in time in the direction opposite to the boost. The number of bubble collisions, 
which is proportional to the spacetime volume in the past light cone of the observer with a 
cutoff at the initial surface, will peak in the direction of the slope. 

It has been shown in [7] that this persistence of memory effect disappears in the limit of 
7 -C 1. The reason is that the number of collisions in this case is dominated by collisions with 
the smallest defects of size ~ Hp^ -C Hp^, whose nucleation centers are in the future of the 
collision time. The same applies to collisions with nucleating defects. This is demonstrated 
by an explicit calculation in Appendix 

VII. DISCUSSION 

We investigated various aspects of collisions between strings and domain walls nucleating 
in an inflating false vacuum and our bubble, which we assume to be expanding in that 
vacuum. These collisions have some similarity to collisions with other bubbles, but in some 
respects they are very different. 

The character of collision between the bubbles crucially depends on whether the vacuum 
energy density of the other bubble is higher or lower than that in our bubble. If it is lower, 
then the invading bubble will expand into our bubble at nearly the speed of light, and all 
objects it encounters on its way, including observers, will be turned into some alien forms 
of matter. If it is higher, then our bubble will expand into the other bubble, and only a 
peripheral region of our bubble will be affected by the collision. 

In the case of defects, the collisions are always rather benign. Defects that do not collapse 
shortly after they enter the bubble, eventually come to rest relative to the comoving observers 
and can be detected through their gravitational effects. We found, however, that there can 
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be no more than a few defects from the multiverse in our observable region, and only if 
their nucleation rate is very high, A ~ 1. Because of their small number, such defects are 
not subject to observational bounds that have been derived for “regular” defects, formed 
in cosmological phase transitions. The tension p of “regular” cosmic strings is bounded by 
[I8] G/r < 10“®. Strings from the multiverse, on the other hand, could have, for example, 
Gp > 10“® without conflict with observational data, provided that the angular extent of 
these strings on the sky is sufficiently small. “Regular” domain walls are essentially ruled 
out by observations, unless their tension is extremely small, cr < (1 MeV)^ [IH]. A wall with 
a larger tension crossing our observable region would introduce a very large gravitational 
perturbation, which would be in conflict with the isotropy of the CMB and of the galaxy 
distribution. The tension of walls from the multiverse can be greater by many orders of 
magnitude, if these walls are located at sufficiently high redshifts. The effect of massive 
defects on the CMB power spectrum has been discussed in [20] for a straight cosmic string 
and in [ 21 ] for a planar domain wall. In particular, it has been argued in [ 21 ] that a massive 
domain wall at about twice the horizon distance from us could account for the hemispherical 
power asymmetry indicated by recent observations. 

Collisions of walls and strings with our bubble can have observational signatures even 
if the defects themselves are not within our present Hubble radius. One can expect the 
observational effects of collisions with domain walls to be similar to the effects of collisions 
with other bubbles. Such collisions will produce round hot or cold spots on the CMB sky and 
a characteristic CMB polarization pattern [SHSlCniE]- Collisions with strings do not have 
rotational symmetry and (unlike collisions with domain walls and with other bubbles) can 
also produce gravitational radiation. We leave a detailed analysis of observational signatures 
of collisions with defects for future research. 

We have estimated the number of potentially observable collisions with defects in our 
past light cone as 

M ~ X{Hf/Ht)\ (94) 

where Hp and Hp are the expansion rates, respectively, of the false vacuum and of the 
slow-roll inflation inside the bubble. This number can be large even if the nucleation rate of 
the defects is A 1. 

We hnally mention some other types of defects which can also nucleate in de Sitter space 
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and collide with our bubble. In a wide class of particle physics models, domain walls can 
be bounded by strings [221 [23] • In such models, the nucleating defects will have the form 
of disc-like domain walls bounded by circular loops of string. The corresponding instanton 
is the usual domain wall instanton with a hole bounded by a string. (The same instanton 
describes nucleation of circular holes in the walls.) 

Monopoles and antimonopoles can be produced in pairs, as discussed in Ref. ns, with 
a separation 2Hp^. Magnetic monopoles behave as localized particles and, considering 
their small numbers, do not produce any signihcant effects. Global monopoles, on the other 
hand, have a non-local distribution of energy outside of their cores and can create strong 
gravitational perturbations m both in our Hubble volume and when they enter the bubble 
in our past light cone. Global monopole and antimonopole in a pair are attracted to one 
another with a force independent of their separation, so pairs that enter our bubble withe a 
separation < will collapse and annihilate. By the same argument as in Sec. IV.G, the 


number of potentially observable collision events is still given by Eq. (94). 
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Appendix A: Coordinate Systems 

In this Appendix we summarize all coordinate systems we used in this paper and the 
relations between them. The false vacuum spacetime can be embedded in a 5-dimensional 
Minkowski space with coordinates (X, lU, U) as 

X^ + W^-V^ = Hf. (Al) 

The surface W = with U > 0 corresponds to a bubble wall with nucleation center at 
X = V = 0. The spacetime of true vacuum inside the bubble can also be embedded in a 5 
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dimensional Minkowski space as 


X2 ^ 1^/2 _y2 ^ 

where W' > and > 0, with the bubble wall at W' = . 

The flat slicing coordinates (t,r, 9,4)) are deflned by 

V = Hp^ sinh Hpt + (A3) 

Xi = 

W = Hp^ cosh HFt - 

with {uz,0Jxi^y) = (cos6^,sin6^cos0,sin6*sin0). And the open slicing coordinates {t,4,9,4>) 
inside the bubble are deflned by 

V = sinh Htt cosh 4 (A4) 

Xi = ^ sinh HtT sinh 

W' = H^^coshHrT, 

In general there is no transformation between {t,r,6,4>) and {t,4,9,4>), since they cover 
different regions of spacetime. However, in the case of Ht = Hp = H, these coordinate 
systems are related in the overlap region through 

coshhfr = coshi^'t — -e^V^ 

2 

sinh TTr sinh .^ = e^V. (AS) 

To represent the whole de Sitter space-time, we use the standard conformal diagram that 
suppresses a 2-sphere at each point and makes use of the conformal version of the closed 
slicing of de Sitter space, namely, 

= tjo ^ 2 X dnl) . (A6) 

H'^sm rj 

The nice thing about this slicing is that it covers the whole space-time, so we can use it to 
represent the evolution of the worldsheet of the defects in different regions by mapping to 
this coordinate system. 

In order to match the space times across the bubble light cone we take —tt < ?7 < 0, 
0 < X < TT for the false vacuum and y — n/2 < rj < arcsin (cosytanh |log 7 |), 0 < y < 7r/2 
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for the true vacuum. The conformal coordinates relate to the flat slicing coordinates that 
cover the false vacuum through 

sin?/ 


r = 


cos T] + cos X 
sinx 


cos 1 ] + cos X ’ 

and relate to the open slicing coordinates that cover the true vacuum through 


(A7) 


sinh HtT = 


sinh {rj + log 7 ) 


, tanh rj = 


smr/ 

cosx’ 


tanh^ = 


smx 

COST] 


(A 8 ) 


Appendix B: Jacobian Matrix 


In this Appendix, we calculate the Jacobian of the transformation matrix from (r„, c) to 
(^m, A0). To simplify the calculation, we dehne 


Then we have 


and 


M = P = cos ( 


^ (1 + M)(1 + P + M(P-1)) 


1 + M2 + P - M2p 


2M 


” (1 + M)(1 + P + M(P-1))’ 

And according to the mapping ( [^ , we have 

2M7 

1 - M2(l - p) + p + 2 MP 7 ’ 

1 + M2(-1 + P) + P + 2 MP 7 


r„. = 


c = 


Thus, 


d{rn,c) 


d (M, P) 

and hnally we obtain 
d (r„, c) 


^(^m, A0) 


(Bl) 


(B2) 


(B3) 


^(1 + M2(-l + P) + pf - 4M2 (-1 + p2) 72 

4M (1 - M2(-1 + P) + P) 72 
((1 + M2(-l + P) + Pf - 4M2 (-1 + P2) 72 )^/^ ’ 

4M (1 - M2(-1 + P) + P) 72 dP 

((1 + M2(-l + P) + Pf - 4M2 (-1 + P2) 72 )^/^ d^m dA(j) 


(B4) 


(B5) 


(B 6 ) 
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Appendix C: Total number of loops within a comoving radius r 


In this appendix, we outline the calculation of the integral (51). We first define new 
variables, A = and C = cAr'^. Then 


dN = dvrA sin amC^rf^drndc, 


= A - + 1 - CfdAdC, 


(Cl) 

(C2) 


and the integration domain is bounded by 


Va 


> 


1, \/A < 1 + VC, c> 0, ^/A>l-y/c. 


(C3) 


One can verify that 0 < \/C < for 1 — \fC < \fA < 1 + \fC^ and < \fC < 
for 1 + ^J~C < ^J~A < VC/r — 1. We also separate the A-integration over 1 + \fC < \/A < 
\fC jr — 1 into two parts, with a separation point at A = 1 + C*. Then the number of loops 
Niv) is expressed as a sum of three integrals. 


»(l+V^)2 


N{r) =Xn^dCl dAWA - ^ (A + 1 - 0)70 

Jo J{i-vcr V 4 

^2 


»(l + x/C)2 I ^ 


4r^ 

‘iA^A-‘-(A+i-cr/c 


■i 

(i+W 


f^+c I i : 

Avr / dC dAJA--{A + l-Cf/C. 

■ / 4 r 2 . V 4 


{l+r)2 


(C4) 

(C5) 

(C6) 


Performing the integrations, we obtain 


N{r) = 


A^2^2(j,2 _|_ 

(l-r2)2 ' 


(C7) 


Appendix D: Collision marks on the sky 

In this Appendix we shall describe the method we used to generate the distribution of 
collision marks on sky. 
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1. Domain walls 


For a domain wall centered on the positive part of the x-axis, the intersection of the 
worldsheets of the wall and the bubble wall can be found from the equations 


— 2 r„r sin 9 cos 0 + = e + e 


-2t 


r = 1 — e 


— t 


(Dl) 


The past light cone of the centered observer intersects the bubble worldsheet at r = ri{To, 7 ), 
where tq is the time elapsed from bubble nucleation to the observation point. Using 


the coordinates dehned in Appendix 


A 


we hnd ri = 5(1 + tanh ^ ) with sinh HtTq = 


sinh(?7o+log7)' 


In the limit of tq S> Hj^^, we have 


1 / , log 7 

r; - 1 — tanh- 

2 V 2 


(D2) 


which reduces to r; 7 ^ for 7 <C 1. The observable part of the intersection is a disc-like 
region bounded by 

1 


7.2 _ g-2t„ 


r = 


2 (r„ sin 9 cos 0 — 1) ’ 


(D3) 


with r < ri, that is 


2(1 — sin 9 cos 0)r; = e 


-2t„ 


-|- 1 . 


(D4) 


If the center of the domain wall is not on the x-axis, the outline of the resulting collision 


mark on the sky can be found from Eq. (D4) by applying a suitable rotation. 

For domain walls with nucleation centers in the blue region of the parameter space, the 
angular extent of this region on the sky is A(j) = 20^, where 


(j)m = arccos 


e‘^^"{2ri — 1 ) — 1 ’ 


(n - l)cOSXn + 7 COS r]n 


CXiL 

ri sin Xn 


(D5) 


For nucleation centers in the green region, the angular extent is A0 = 27r — 20^- “Blue” 
domain walls always have A0 < tt. This means that the area affected by the collision is less 
than half of the sky. For green walls, A0 can either be greater or smaller than vr, but we 
hnd numerically that collisions with A 0 > tt are much more likely. 

An example of the parameter space (c = e*”, r„) for r; = 0.6, = 2 is shown in Fig. 

As before, nucleation centers in the blue area correspond to walls hitting the bubble from 


16 
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FIG. 16: The parameter space of domain walls with n = 0.6 and Se = 2. Both green and blue 
regions continue to infinity in the direction of r„. 


the right, and nucleation centers in green areas (both dark and light) correspond to walls 
hitting the bubble from the left. The dark and light green regions correspond to collisions 
with A0 > TT and A0 < vr, respectively. The curve B is the bubble boundary. 

For a collisions to be observable, we require that tc < ti, which gives 


1 

^/{rn - l)(r„ - 2ri + 1) 

for blue walls and 


(D6) 


1 

V'(r„ + l)(r„ + 2ri - 1) 

for green walls. The resulting cutoff curves are labeled C and D 
We also require tf < tc for blue walls and tf < 0 for green 
constraints are, respectively. 


(D7) 

in the hgure. 

walls. The corresponding 


rn> I- 




(D8) 


c 

and c < Se- They exclude the gray region of the parameter space in the hgure. 

The blue and green areas of the parameter space extend to inhnite values of r„, but the 
volume of the spacetime region that they specify is hnite. We randomly pick parameter sets 
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{tn,'f'n) from the allowed region according to the distribution (42) 


dN = AnXr'^drnC^dc. (D9) 

For each parameter set we choose a random direction and plot the resulting collision mark 
on the sky with the Mollweide projection. 


2. Cosmic Strings 


The intersection between the worldsheets of a string loop centered on the x-axis and the 
bubble wall can be found from the equations 


+ 1 


2(1 — r„sin6*cos0)r = e 
(sin 6 cos 0 tan a — cos 6)r = Vn tan a 


(DIO) 


The projection of the intersection line on the celestial sphere can be easily expressed in a 
parametric form ^(r), 0(r); we do not reproduce the explicit expressions here. An observer 
at r = 0 can see only the part of this line at r < r;. If the center of the loop is not on the 
x-axis, the resulting collision mark on the sky can be found by applying a suitable rotation. 

As before, we impose the conditions tj < < U for all loops and tf < 0 for loops 


with nucleation centers in the green region. We then use the distribution (49) to randomly 


choose the nucleation parameters {tn,rn,a) in this restricted parameter space. Finally, we 
randomly pick a direction to the loop nucleation center and plot the collision mark with the 
Mollweide projection. 


Appendix E: Distribution of collision events for a non-central observer 

We shall now study collision events within the past light cone of a non-central observer. 
It will be convenient to introduce flat-slicing coordinates inside the bubble, 

= dt'"^ — [dx'^ + + dz'‘^) , (El) 

in terms of which the past light cone of an observer at (to,XQ, 0, 0) is 

laf — XqI = = -sj {x' — (E2) 
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and the bubble wall is 


If' 


= 1 — e = \Jx'"^ + y’’^ + 


,12 


(E3) 


The intersection of the past light cone and the bubble wall can be found by combining the 
above equations, 

1 Xq + 1 — — 2x'q cos -0 


1 — e ^^*0 — Xq cos -0 


1 (1 — e 


(E4) 

—Hrt' / T’ 

/ 1 — e — Xq cos Ip 

where we have introduced cos'p = x'/r' and r' = |f'|. In the limit of I'q —)■ oo, the spatial 
projection of the intersection surface is an ellipsoid 

which has one focus at the origin and the other at the observer’s location Xq. Note that here 
Tg is a comoving coordinate. 

In order to hnd the number of string loops that intersect with this surface, we should hrst 
rewrite the ellipsoid in terms of the flat-slicing coordinates outside the bubble. The junction 
between the interior and exterior spacetimes requires that the physical radii of two-spheres 
match at the location of the wall, namely 




r . 


For points on the bubble wall, this gives 


Te. = 


(E7) 


(E8) 


(l-7)r; + 7' 

According to Ref. ca, nucleating string loops are uniformly distributed on a constant t 
surface. The number of loops of physical radius R whose centers are located in a 3-volume 
element dV = 27rrpj^y sin 'ipdrp^yd'ip is 


dN = \^dV. 


(E9) 


This distribution is cut off at R ~ Hp^. 


In order to intersect with the ellipsoid, a loops of radius R <C Hp ^ should have its center 
in the range 


H, 


-1 


- -F Hp^ < Tphy < H, 


1 — r 


-1 


1 — r 


j + R^ 


(ElO) 
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where we used e = 1 — r'^. Then the distribution of intersections with respect to ^|J is 


dM 

d'lp 


~ 2'kX sinh’ / —- 

^ ' -1 


'HZ 




_j ^phy'^^phy 


^ X -2 • , Re 

~ — A 7 sm'^ 


/2 


= Ja7 


-2 


(1 - r ^)2 

(xq^ — 1)^ sinV’ 


3 ' (1 + Xq — 2x'q cos ' 

where we have evaluated the integrals in the limit of 7 1 . 


(Ell) 

(E12) 

(E13) 


Introducing the zenith angle 'ipo with respect to the non-central observer, the above 
distribution becomes 


dAf 

di>o 


~ 27rA7 ^ sin -00 


(E14) 


We thus see that there is no anisotropy in the number distribution of collision events in the 
limit of 7 -C 1. The total number of collision events within the past light cone can be found 


by integrating Eq. (E14) over 


V = yA7-^ 


(El 5) 


which is the same as (82). 
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